Nodal structure of quasi-2D superconductors probed by magnetic field 
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We consider a quasi two-dimensional superconductor with line nodes in an in-plane magnetic 
field, and compute the dependence of the specific heat, C, and the in-plane heat conductivity, k, 
on the angle between the field and the nodal direction in the vortex state. We use a variation of 
the microscopic Brandt-Pesch-Tewordt method that accounts for the scattering of quasiparticles off 
vortices, and analyze the signature of the nodes in C and k. At low to moderate fields the specific 
heat anisotropy changes sign with increasing temperature. Comparison with measurements of C 
and k in CeCoIns resolves the contradiction between the two in favor of the d x 2_ y 2 gap. 

PACS numbers: 74.25.Fy, 74.20.Rp, 74.25.Bt 
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Introduction. Remarkable properties of superconduc- 
tors below the transition temperature, T c , are due to for- 
mation of the Cooper pairs and opening of the energy 
gap in the single particle spectrum. The structure of the 
gap is intimately related to the nature of the pairing in- 
teraction, and therefore its experimental determination 
is of much importance. In simple metals the phonon- 
mediated electron attraction and the gap are isotropic 
around the Fermi surface. In contrast, in the majority of 
recently discovered superconductors the gap is strongly 
anisotropic, and often vanishes (has nodes) for selected 
directions in the momentum space. 

While the existence and topological structure (line vs. 
point) of the nodes can be inferred from the power laws in 
the temperature dependence of the thermodynamic and 
transport properties, experimental determination of the 
nodal directions is more challenging. Currently the most 
widely used probes utilize an applied magnetic field, H, 
as a directional probe of unpaired quasiparticles (QP). 

Low energy quasiparticles only exist in the near-nodal 
regions. They are excited by temperature, T, at all 
nodes, and by the field, predominantly at the nodes 
where the QP velocity is normal to H . Consequently the 
QP density changes with the relative orientation of the 
field with respect to the nodes. It was predicted that the 
low-energy density of states (DOS), and the electronic 
specific heat, C, have minima when the magnetic field is 
aligned with the nodes 0. 

Due to challenges in measuring the electronic contribu- 
tion to C(T, H), so far experimental results exist for few 
materials 0, Q ■ At the same time the anisotropy of the 
electronic thermal conductivity, under a rotat- 

ing field was determined in several systems 0, 0, @ ■ The 
oricntational dependence of the in-plane k(H) is more 
complex than that of the specific heat as it combines 
the anisotropy due to nodal structure with that due to 
the difference in scattering normal to and parallel to the 
vortices. Theoretical interpretation of the anisotropy of 
transport properties proved elusive as the existing theo- 
ries have difficulties accounting for the two modulations, 
and cannot unequivocally determine whether a local min- 



imum or a maximum occurs when the field is along a 
node. Consequently, more theoretical work is needed to 
examine the experimentally conjectured gap anisotropy 
00. Moreover, in heavy fermion CeCoIns, where data 
on both C and k anisotropy exist, the conclusions appear 
contradictory: d xy (d x 2_ y 2) gap symmetry was inferred 
from C (k) Hi]- 

In most approaches the effect of H on the QPs is in- 
cluded via the Doppler energy shift due to supercurrents 
around the superconducting vortices , and the QP life- 
time is not affected directly 0- In reality, vortices scat- 
ter quasiparticles carrying the heat current and not just 
shift their energy. Consequently, the behavior of k(T, H) 
is determined by the competition between the enhance- 
ment of the DOS and the vortex scattering. There are in- 
dications that vortex scattering affects transport already 
at moderate, compared to the upper critical field H C 2, 
fields 0, , but its effect for the in-plane field and on 
the anisotropy in C(H) and k(H) has not been studied. 

In this Letter we present a unified microscopic ap- 
proach to computing the anisotropy of thermodynamic 
and transport properties of nodal superconductors in the 
vortex state at moderate to high magnetic field. We ap- 
ply the method to quasi two-dimensional superconduc- 
tors with line nodes, and consider the field rotated in the 
xy plane. We are able to account for both twofold and 
fourfold variations of n(T,H) with angle. For d x z_ y 2 
gap, the competition between the transport scattering 
rate and the DOS leads to switching from minima to 
maxima for H along the nodes in the fourfold part of 
k(T, H) upon increasing temperature. We find that, due 
to the field dependence of the single particle lifetime, the 
minima and the maxima in the specific heat also switch 
at higher T and H. Hence in a wide T-H range the max- 
ima (rather than minima) in C(T,H) indicate a nodal 
direction. Our results for the d x 2_ y 2 gap are in a semi- 
quantitative agreement with experiment on CeCoIns. 

Model and approach. We consider a quasi-two dimen- 
sional system with a model Fermi surface (FS) p 2 = 
Px + P y ~ (f 2 P 2 f)cos(2sp z /r 2 pf). The parameters, r 
and s, determine the corrugation amplitude along the 
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z-axis and the ratio Vf lX /vf t ± (at p 2 + p 2 = p 2 ) re- 
spectively. We assume a separable pairing interaction, 
y(Pi p') = Voy(p)y(p'), where y(p) are the normalized 
basis functions for the angular momentum eigenstates, so 
that y(p) = v / 2cos2</>p (y(p) = v / 2sin20p) for d x i^ y i 
(d xy ) gap. Here <fip is the angle between the projection of 
the momentum p onto the xy plane and the x-axis. The 
magnetic field, H, is in the xy plane, at an angle <po to 
the a;-axis. We ignore Zeeman splitting, which is justified 
for most systems, see below. 

We rotate x and y axes around z to choose the 
field direction as the new x-axis, and model the spa- 
tial dependence of superconducting order parameter 
by an Abrikosov-like solution, A(R, p) = A(R)3^(p), 
where A(R) = £„ A »( R l n >- Here ( R l n > = 
XX f n ((z - A 2 fc,)/ViA) exp(z/c y y), $ n (z) is the 
eigenfunction of the n-th state of a harmonic oscillator, 
A 2 = hc/2\e\H, and normalized cjg {J2 ky \C^\ 2 = 1) 
determine the structure of the vortex lattice. The factor 
■s/s ensures that the functions <!>„ are the solutions to the 
linearized equations for the order parameter |lllll2j . and 
A n is the amplitude of the nth spatial component. In s- 
wave superconductors the vortex state is a superposition 
of n = states only; however, for singlet pairing in the 
higher angular momentum channels the self-consistency 
requires admixture of other even n |l3l flij . In our cal- 
culations we use n = 0,2,4, since inclusion of higher n 
does not alter the results. 

We use the quasiclassical Green's function approach 
[l5| . where the equation for the retarded (index R) 
anomalous (Gor'kov) Green's function, is given by 



-2£e + v/(p) V 



/ fl (p,R;e) 

2Aig R (p,R;e). (1) 



Heree = e-T, R (R;s), A = A(p,R) + A^ p , and we treat 
the self-energies A R mp due to impurity scattering 
in the self-consistent T-matrix approximation, and focus 
here on the unitarity limit (scattering phase shift tt/2) 
for clean systems (normal state scattering rate T <C T c ). 
The normal electron Green's function g R (p, R; e) is de- 
termined via the normalization condition {g R ) 2 — f R f R = 
-7T 2 , where / fl (p,R;e) = f R (-p, R; -e)*. Eq.© is 
complemented by the self-consistency condition for A. 

To solve the quasiclassical equations we make use of 
the approximation due to Brandt, Pesch, and Tewordt 
(BPT) 0, and replace g R (R,p,e), by its spatial aver- 
age. If K is a vector of the reciprocal vortex lattice, the 
Fourier components g R (K) cx exp(— A 2 K 2 ), hence the 
spatial average K = dominates [lfj. The method is 
nearly exact at H < H C 2, but gives semi-quantitatively 
correct results down to much lower fields. In extreme 
type-II 01 j an d in the nodal superconductors 0, 0] 
the method remains valid over almost the entire range 



H c i <C H < H c2 0, and was used to study unconven- 
tional superconductors in the vortex state ^ [ill [2(| • 

With averaged g R we solve Eq.Q by expanding 
A,f R ,f R in the orthonormal set {(R|n)}, and using 
the ladder operators for the oscillator states, a and , 
[Hill to rewrite v/(p) (V R - i'(2|e|/ftc)A(R)) = u_a- 
v+a\ where v± = [v z j ^/s±iv y ^/s\/{A^/2). Lengthy but 
straightforward calculation gives the closed form expres- 
sions for functions f R , f R , and g R for a set A„, which is 
then determined sclf-consistcntly. For n = only we find 



9 R (p;e) 



(2) 



'y 2 (p)w(^) 



similar to Refs.0,0. Here W(z) = exp(— z 2 ) erfc(— iz), 
and the dependence on the field direction is via the 
rescaled component of the Fermi velocity normal to H , 
\vf (P/)| 2 = v fz/ s + sv fy ^ ne expression for arbitrary 
n has a similar structure, but involves combinations of 
A„ and higher derivatives of W(z) (l2^ . 

The density of states (DOS) is N(T,H;e)/N f = 
( — 1/n) J dpf 3m g R {p; e), where Nf is the normal state 
DOS per spin at the Fermi level. From the DOS we deter- 
mined the entropy, S(T, H), computed the specific heat, 
C = TdS/dT, by numerical differentiation, and verified 
that far from the transition T C (H) the expression 



C(T,H) = 



ds 



s 2 N(T, H; e) 
T 2 cosh 2 (e/2T) 



(3) 



gives accurate results. We use Eq.© below. 

The thermal conductivity tensor can be obtained us- 
ing Keldysh technique [l2T | , and has a particularly simple 
form for the Born and unitarity scattering limits |22l |23| 



« -(T,g) 
T 



de e 2 , n £ 

7T cosh — 

2T T 2 2T 



(4) 



x J dpv f>i v f!j N(T,H;p,s) T H (T,H;p,£) . 

Here N is the angle-dependent DOS, and th has the 
meaning of the transport lifetime due to both impurity 
and vortex scattering. For n = channel [9l.llll.l24j. 
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combination of the IF-function and its derivatives [12j . 

Results. We show the results for s = 0.5, r = 0.5, 
which gives ifff? /H% re 2.45 similar to CeCoIns, and for 
T/{2ttT c ) = 0.007 Q. The main qualitative difference 
from previous results is already clear from the behavior of 
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FIG. 1: (a) Low energy DOS for H\\ antinode and H \\node at 
T = 0.05 T c . Arrows indicate e(H), see text, (b)-(d) Angle- 
resolved DOS plotted at each position on the FS. Net DOS is 
proportional to the red shaded area, see text. 



the DOS shown in Fig^a). In the Doppler shift method 
the DOS anisotropy, N(e, H\\antinode) — N(s, H\\node), 
increases as s/H at e = 0, and vanishes at e ~ VFy/s/K 
Pl|2(|. In contrast, the anisotropy in the residual (e = 0) 
DOS has a maximum at H ~ 0.1-ff C 2, and reverses above 
the field H ~ Q.5H C 2- Below this field the anisotropy also 
changes sign at a finite e(H), see Fig^a). 

Since W(z) and its derivatives in Eq.(J2J) are complex 
functions, our N(e) cannot be obtained from the BCS 
DOS by a simple energy shift: there is an anisotropic sin- 
gle particle scattering rate due to scattering from vortices 
(vanishingly small along H, largest normal to H). This 
occurs since in the BPT method g R is averaged incoher- 
ently in different unit cells of the vortex lattice. Vortex 
scattering is pairbreaking, and hence enhances the DOS. 

Consider the contribution to the DOS from different 
FS regions and focus first on e = 0, Fig^b),(c). At low 
fields the vortex scattering rate is small, and the unpaired 
states emerge only near the nodes. When H\\node the 
number of such states at the nodes aligned with the field 
is small, while H\\antinode produces states at all nodes, 
see FigHJc), and N(0, H\\ antinode) > N(0,H\\node) as 
in the Doppler shift method Q . At higher fields the pair- 
breaking is stronger, and, for H \\antinode, it generates 
the unpaired states in all directions on the FS except close 
to the field direction, see Fig^b). For if || node the field- 
induced states appear everywhere on the FS [ Fig^b)], 
leading to the anisotropy reversal • 

Now consider < e -C A , Fig^d). In a pure system 
at zero field the dominant contribution to N{e) is from 
sharp peaks at directions p e , close to the nodes, such 
that |A(p e )| = e. Scattering redistributes the spectral 
weight from the peak. Vortex scattering at p 6 , increases 
with H, and is stronger for H _L p e than for JJ||p e . For 
e > e(H) H\\node fills the near- nodal states, but does 



FIG. 2: Specific heat anisotropy for different T at two fields. 
Nodal direction cj>o = 45°. Left panel: H — 0.26H C 2, inver- 
sion of the anisotropy occurs at To ~ 0.15T C . Right panel: 
H = 0.527T C 2, no minimum for the field along the node at 
any T. Some curves are shifted by the amount indicated in 
parentheses. 



not broaden the peak nearest to p e significantly, Fig[IJd). 
For H\\antinode most of the weight in the peak is shifted 
away. This leads to a higher DOS at e > e(H) for the 
field along the node, see Fig^a). 

Fig. |21 shows the specific heat anisotropy. The domi- 
nant contribution to C(T, H) is from the DOS at energies 
e ~ 2.4T, see Eq.©. Therefore even at low fields, when 
iV(0, if |j antinode) > N(0, H \\node), the anisotropy of 
the specific heat is reversed at T Q - 0.1-0.2 T c 0. Be- 
low To, C(T, H) has a minimum when the field is applied 
along a nodal direction, while at higher T it has a maxi- 
mum for this field orientation. At higher fields H along 
a node gives a maximum in the angle-dependent specific 
heat at all but the lowest T. This result directly affects 
the experimental determination of the nodal directions. 

The measurements on CcCoIns were carried out for 
0.18 < H/H c2 < 0.5, and for T > 0.1T C 0. For 
this system H C 2 is Pauli limited [2^|, and the orbital 
critical field, , may be as high as 2.bH C 2- We 



find that Tq{H) in Fig. 2 is weakly field dependent for 
0.1 < H/H^ rh) < 0.3. Therefore our results indicate 
that maxima, rather than minima in C (T, H) occur for 
H\\node in this regime, and the data of Ref. support the 
d x 2_ y 2 gap symmetry (rather than the d xy order inferred 
by the authors from the low T, low H theory 0). While 
the BPT approach likely overestimates the vortex scat- 
tering at low H/H C 2, the extended range of this shape of 
C(T,H) is also in favor of d x 2_ y 2 pairing. 

This conclusion is supported by the analysis of the 
thermal conductivity, shown in Fig. for the heat cur- 
rent Comparison can be made with k(T,H) as 
shown in Ref. for T < 0.5T C and fields H < 0.3H c2 . 
In this range we find the overall shape of the curve and 
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FIG. 3: Thermal conductivity anisotropy for different T at a 
H = 0.26H c2 (left panel), and for different H at T = 0.25T C 
(right panel). Curves are shifted for clarity by the amount 
indicated. Left panel: notice the reversal of the 0° — 90° 
anisotropy and the change from minimum to maximum at 
45° . Right panel provides direct comparison with Ref|^. 



the amplitude of the peak (~ 1-3%) at the nodal angle 
always similar to those in the right panel Fig. [3J We also 
confirmed that d xy gap is inconsistent with the results of 
Ref. 0. The orientational dependence of C(T, if) for the 
d xy gap is obtained by a 45° degree rotation in Fig[21 In 
contrast, the variation of k(T, H) with angle is different 
for the two cases since the heat current j q \\x is alon g th e 
antinodal (nodal) direction for the d x i_ y i (d xy ) gap |12j . 

The fourfold angle dependence of n{H) is superim- 
posed on the twofold variation due to relative orienta- 
tion of j q and H. The field is strongly pair breaking for 
p _L H, and QPs which contribute the most to k(T, if), 
are created for H _L j q [3(|. Consequently, at low T, 
k(T,H _L j q ) > n{T,H\\j q ). At higher T the QP arc 
thermally excited, and k(T,H _L j q ) < n(T, H\\j q ) due 
to vortex scattering, see Fig. k(T, if \\node) has a lo- 
cal minimum at low T, H, but develops a maximum at 
higher T. At yet higher T, k(T, if) is essentially twofold. 

Conclusions. We developed a fully self-consistent 
microscopic calculation of the magnetic field induced 
anisotropy of thermal and transport properties of nodal 
superconductors. We applied the method to quasi-two 
dimensional systems with line nodes. We 1) found that, 
while at very low temperatures and fields the specific 
heat anisotropy is in agreement with the Dopplcr shift 
analysis, at higher T and H it changes sign, and, over a 
wide range of fields and temperatures, exhibits a maxi- 
mum, rather than a minimum, for the field applied along 
a nodal direction; 2) accounted simultaneously for the 
two-fold and fourfold pattern in the angle-dependence of 
the thermal conductivity; 3) used the approach to re- 
solve a controversy regarding the symmetry of the order 
parameter in CeCoIns in favor of d x 2_ y 2 gap. We believe 



that our approach provides a reliable tool to be used 
alongside experimental measurements for determination 
of the nodal directions in novel superconductors. 
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